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Frequency shifts of the 199 Hg+ 5d 10 6s 2 S 1/2 (F = 0, M F = 0) to 5d 9 6s 2 2 D 5/2 (F = 2, M F = 0) 
electric-quadrupole transition at 282 nm due to external fields are calculated, based on a com- 
bination of measured atomic parameters and ab initio calculations. This transition is under 
investigation as an optical frequency standard. The perturbations calculated are the quadratic 
Zeeman shift, the scalar and tensor quadratic Stark shifts, and the interaction between an ex- 
ternal electric field gradient and the atomic quadrupole moment. The quadrupole shift is likely 
to be the most difficult to evaluate in a frequency standard and may have a magnitude of about 
1 Hz for a single ion in a Paul trap. 

Key words: atomic polarizabilities; electric quadrupole interaction; mercury ion; optical fre- 
quency standards; Stark shift; Zeeman shift. 



1. INTRODUCTION 

It has long been recognized that a frequency stan- 
dard could be based on the 282 nm transition between 
the ground 5d 10 6s 2 Si/2 level and the metastable 5d 9 6s 2 
2 D 5 / 2 level of Hg + [[jj. The lifetime of the upper level 
is 86(3) ms Q, so the ratio of the natural linewidth Au 
to the transition frequency vq is 2 x 10~ 15 . (Unless oth- 
erwise noted, all uncertainties given in this paper are 
standard uncertainties, i.e., one standard deviation esti- 
mates.) Doppler broadening can be avoided if the tran- 
sition is excited with two counter-propagating photons, 
as originally proposed by Bender et al. O and subse- 
quently demonstrated by Bergquist et al. 151. However, 
optical Stark shifts are greatly reduced if the transition 
is driven instead with a single photon by the electric- 
quadrupole interaction. In this case, Doppler broaden- 
ing can be eliminated if the ion is confined to dimen- 
sions much less than the optical wavelength, as was first 
demonstrated by Bergquist et al. Q. 

Recently, the (F = 0, M F = 0) to (F = 2, M F = 0) hy- 
perfine component of the 199 Hg + 5d 10 6s 2 Si/ 2 to 5d 9 6s 2 
2 D5/2 single-photon transition has been observed with a 
linewidth of only 6.7 Hz by Rafac et al. ||. A laser servo- 
locked to this transition is an extremely stable and repro- 
ducible frequency reference. New developments in opti- 
cal frequency metrology ||, |?| may soon make this system 
practical as an atomic frequency standard or clock. 

While the (F = 0, M F = 0) to (F = 2, M F = 0) hy- 
pcrfinc component has no linear Zeeman shift, it does 
have a quadratic Zeeman shift that must be accounted 
for. In addition, there is a second-order Stark shift and a 
shift due to the interaction between the electric-field gra- 
dient and the atomic electric-quadrupole moment. None 
of these shifts has yet been measured accurately, so it is 
useful to have calculated values, even if they are not very 



precise. Also, it is useful to know the functional form of 
the perturbation, even if the magnitude is uncertain. For 
example, the quadrupole shift can be eliminated by av- 
eraging the transition frequency over three mutually or- 
thogonal magnetic-field orientations, independent of the 
orientation of the electric-field gradient. 



2. METHODS AND NOTATION 

The quadratic Zeeman shift can be calculated if the hy- 
pcrfme constants and electronic and nuclear ^-factors are 
known. Similarly, the quadratic Stark effect can be cal- 
culated from a knowledge of the electric-dipole oscillator 
strengths. The quadrupole shift depends on the atomic 
wavefunctions. Some of these parameters have been mea- 
sured, such as the hyperfine constants and some of the 
oscillator strengths. There are also published calcula- 
tions for some of the oscillator strengths. 

Here, we estimate, by the use of the Cowan atomic- 
structure codes, values for parameters for which there 
are neither measured values nor published calculations. 
The Cowan codes are based on the Hartree-Fock ap- 
proximation with some relativistic corrections ||. The 
odd-parity configurations included in the calculation 
were 5d 10 np (n = 6,7,8,9), 5d 10 5/, 5d 9 6s6p, 5d 9 6s7p, 
5d 9 6s5/, and 5d 8 6s 2 6p. The even-parity configurations 
were 5d 10 ns (n = 6,7,8,9,10), 5d w nd (n = 6,7,8,9), 
5(i 9 6s 2 , 5d 9 6s7s, 5d 9 6s6d, and 5e? 9 6p 2 . Recently, San- 
sonetti and Reader have made new measurements of the 
spectrum of Hg + and classified many new lines |j] . They 
also carried out a least-squares adjustment of the energy 
parameters that enter the Cowan-code calculations in or- 
der to match the observed energy levels. We use these 
adjusted parameters in our Cowan-code calculations. 

As one test of this method of calculation, we estimated 
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the weakly allowed 10.7 //m 5d w 6p 2 P 1/2 to 5d 9 6s 2 2 D 3/2 
electric-dipole decay rate. This decay is allowed only be- 
cause of configuration mixing, since it requires two elec- 
trons to change orbitals. The calculation shows the de- 
cay to be due mostly to mixing between the 5d 10 6p and 
5d 9 6s6p configurations. The calculated rate is 111 s _1 ; 
the measured rate is 52(16) s _1 Q. Another test is the 
electric-quadrupole decay rate of the 5d 9 6s 2 2 D 5 / 2 level 
to the ground level. The calculated rate is 12.6 s" 1 , and 
the measured rate is 11.6(0.4) s _1 . Similar calculations 
have been carried out by Wilson [|l0| . 

Let Hq be the atomic Hamiltonian, exclusive of the hy- 
perfine and external field effects, which are treated as per- 
turbations. For convenience, we denote the eigenstates 
of Hq corresponding to the electronic levels 5d 10 6s 2 Si/ 2 
and 5d 9 6s 2 2 D 5 / 2 having J z eigenvalue Mj by |S 1/2 Mj) 
and |D 5/2 Mj), respectively. 

The corresponding eigenvalues of Hq are denoted 
W(S,l/2) and W(D,5/2). An arbitrary eigenstate of 
Hq with eigenvalue IF (7, J) and electronic angular mo- 
mentum J is denoted I7 J Mj). Since 199 Hg + has in 
addition a nuclear angular momentum J, where / = 1/2, 
the complete state designation is \jJFMp), where F is 
the total angular momentum, and Mp is the eigenvalue 
ofF z . 

3. QUADRATIC ZEEMAN SHIFT 

In order to calculate the energy shifts due to the hy- 
pcrfine interaction and to an external magnetic field 
B = Bz, we define effective Hamiltonian operators Hg 
and H D that operate within the subspaces of hyper- 
fine sublevels associated with the electronic levels 5d 10 6s 
2 S]y 2 and 5d 9 6s 2 2 D 5 / 2 respectively: 

= hA s I ■ J + 3j(S)/x b J ■ B + g'jp B I ■ B, (1) 

H D = hA B I ■ J + g.j(T>)p B J ■ B + g' lt i B I ■ B, (2) 

where Ag and Ad are the dipole hyperfine constants, 
gj(S) and <?j(D) are the electronic g- factors, g\ is the 
nuclear g-factor, h is the Planck constant, and is 
the Bohr magneton. All of the parameters entering Hg 
and H D are known from experiments, although a more 
accurate measurement of gj(D) would be useful. The 
ground-state hyperfine constant As has been measured in 
a 199 Hg + microwave frequency standard to be 40 507.347 
996 841 59 (43) MHz fnj. The excited-state hyperfine 
constant Ad has been measured recently by an extension 
to the work described in Ref. ||, in which the differ- 
ence in the frequencies of the |S 1/2 0) to |D 5/2 2 0) 
and the \S 1/2 0) to |D 5/2 3 0) transition frequen- 
cies was determined to be 3A D =2 958.57(12) MHz |J|, 



in good agreement with an earlier, le ss p recise measure- 
ment by Fabry-Perot spectroscopy [pL3|| . The ground- 
state electronic <?-factor gj (S) was measured in 198 Hg+ 
by rf-optical double resonance to be 2.003 174 5(74) jjj. 
The excited-state electronic g-factor <?j(D) was measured 
in 198 Hg + by conventional grating spectroscopy of the 
398 nm 5d 10 6p 2 P 3/2 to 5d 9 6s 2 2 D 5/2 line to be 1.198 
0(7) Jill- The difference in gj(S) or gj(D) between 
i98jjg+ anc j i99jjg+ j g estimated to be much less than 
the experimental uncertainties. The nuclear g-factor g'j 
is -5.422 967(9) x 10~ 4 0. The measurement was made 
with neutral ground-state 199 Hg atoms, so the diamag- 
netic shielding factor will be slightly different from that 
in the ion. However, this is effect is negligible, since the 
magnitude of g\ is so small compared to <?j(S) or gj (D). 

The determination of gj (D) could be improved by mea- 
suring the optical-frequency difference between two com- 
ponents of the 282 nm line and the frequency of a ground- 
state microwave transition at the same magnetic field. 
Since the uncertainty in the quadratic Zeeman shift is 
due mainly to the uncertainty in <?j(D), it is useful to see 
how accurately it can be estimated theoretically. The 
Lande (/-factor for a 2 D 5 / 2 state, including the correction 
for the anomalous magnetic moment of the electron, is 
1.200 464. The Cowan-code calculation shows that the 
configuration mixing does not change this value by more 
than about 10~ 6 , i.e., 1 in the last place. There are sev- 
eral relativistic and diamagnetic corrections that modify 
gj(D), one of which, called the Breit-Margenau correc- 
tion by Abragam and Van Vleck ]l7[| , is proportional to 
the electron mean kinetic energy. The other corrections 
are more difficult to calculate. The Cowan-code result for 
the mean kinetic energy of an electron in the 5d orbital of 
the 5d 9 6s 2 configuration is T = 19.32 hcR^o, where i?oo 
is the Rydberg constant. Using this value, we obtain a 
theoretical value of g.j(D), including the Breit-Margenau 
correction, of 1.199 85, which disagrees with the the ex- 
perimental value by 1.85 x 10~ 3 , which is 2.6 times the 
estimated experimental uncertainty of Ref. fl5| . If we 
calculate gj(D) for neutral gold, which is isoelectronic 
to Hg + , by the same method, we obtain a value which 
differs from the accurately measured experimental one 
H by (7 ± 2) x 10~ 5 . Thus, the error in the calculated 
value for gj (D) of 199 Hg + might be less than 1 x 10~ 4 , 
but it is impossible to be certain of this, since there are 
uncalculated terms. Measurements of the 199 Hg + optical 
clock frequency at different values of the magnetic field 
should result in a better experimental value for gj(D) in 
the near future. 

For low magnetic fields (B less than 1 mT), it is suffi- 
cient to calculate the energy levels to second order in B. 
To this order in B, the energies of the hyperfine- Zeeman 
sublevels for the ground electronic level are 
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3hA s [ 9J (S)-^] 2 M |S 2 



W(S, 1/2, 0,0,3) = W(S, 1/2) - 
W(S, 1/2, 1,0,3) = W(S, 1/2) + 
W r (S,l/2,l,±l,B) = W(S,l/2) + 



4 4/iA s 

[g 7 (S) - g'j] 2 >jg 
4 4/iAs 
^4g [gj(S) +g^ B 3 



For the 5ei 9 6s 2 2 D 5 / 2 level we have 



W(D,5/2,2,0,B) 


= W(D,5/2)- 


W(D,5/2,2,±1,£) 


= W(D,5/2)- 


W(D,5/2,2,±2,£) 


= W(D,5/2)- 


PV(D, 5/2, 3,0,B) 


= W(D,5/2) + 


W(D,5/2,3,±1,E) 


= W(D,5/2) + 


W(D,5/2,3,±2,S) 


= PF(D,5/2) + 


W(D,5/2,3,±3,S) 


= W(D,5/2) + 



Here, W(7, J, F, Mp,B) denotes the energy of the state 
|7 JFMp), including the effects of the hyperfine interac- 
tion and the magnetic field. 

At a value of B of 0.1 mT, the quadratic shift of the 
|S 1/2 0) to |D 5/2 2 0) transition (optical clock tran- 
sition^ — 189.25(28) Hz, where the uncertainty stems 
mainly from the uncertainty in the experimental value 
of gj (D). In practice, the error may be less than this if 
the magnetic field is determined from the Zeeman split- 
tings within the |D 5/2 F Mp) sublevels. The reason 
is that an error in <?j(D) leads to an error in the value 
of B inferred from the Zeeman splittings, which partly 
compensates for the <?j(D) error. If instead we use the 
calculated value of g.j(D), the quadratic shift for B = 0.1 
mT is —189.98 Hz, where the uncertainty is difficult to 
estimate. 



4. QUADRATIC STARK SHIFT 

The theory of the quadratic Stark shift in free atoms 
has been described in detail by Angel and Sandars [ jl9| . 
The Stark Hamiltonian is 



(3) 
(4) 
(5) 

(6) 
(7) 
(8) 
(9) 
(10) 

(11) 

where fi is the electric-dipole moment operator, 

/Lt = -e^V;, (13) 

i 

and E is the applied external electric field. In Eq. (|l3|), 
r*i is the position operator of the ith electron, measured 
relative to the nucleus, and the summation is over all 
electrons. 

First consider an atom with zero nuclear spin, such as 
i98jjg+ rp Q geconc [ orc [ er m the electric field, the Stark 
shifts of the set of sublevels \jJMj) depend on two pa- 
rameters, a sca i ar (7, J) and extensor (7 > J)t called the scalar 
and tensor polarizabilities. In principle, when both mag- 
netic and electric fields are present but are not parallel, 
the energy levels are obtained by simultaneously diago- 
nalizing the hyperfine, Zeeman, and Stark Hamiltonians. 
In practice, the Zeeman shifts are normally much larger 
than the Stark shifts, so that Hp does not affect the di- 
agonalization. In that case, the energy shift of the state 
\^JMj) due to Hp is 
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H E = -fi- E, 



AW E (l,J,Mj,E) 



(12) 

-5 a scalar(7 5 J)E 2 - 



rate 



[3M 2 - J(J+1)] 
J(2J- 1) 



(3F 2 - E 2 



(14) 
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Treating Hp by second-order perturbation theory leads to the following expressions for the polarizabilities |19 

^scalar 

ft, J) 



8 7 re _ V \(jJ\\^\h'J')\ 2 



3(2J+l)^M/( 7 ',J')-W / (7^) : 



ate 



c (7, J) = 87T60 



10J(2J-1) 



3(2J + 3)(J + 1)(2J+ 1) 



1/2 



D-D-' 

7',/' k ' 



\hJ\\^\h'Jl)f 



J J J' J W(7', J') - W(7> J) ' 



(15) 
(16) 



The summations are over all levels other than \jJ). Equations ( |l5|) and(16) can be rewritten in terms of the oscillator 
strengths / 7 ,/. 7 'j': 



( 47re e 2 ^ 2 ^ 

a S calar(7i J) = ~ / , 



f- 



J.J,f'J' 



me fa [T^( 7 ',J')-^( 7 ,J)] 2 ' 



4-KCoe h 



2fc2 r 



30J(2J- 1)(2J + 1) 



(2J + 3)(J+1) 



J 7 /J/ v > 



ftJn'J' 



j j j' I [w(y,j')-w(^,JW 



(17) 
(18) 



where m e is the electron mass. The tensor polarizability is zero for levels with J < 1, such as the Hg + 5d 10 6s 2 Si/ 2 
level. 

For an atom with nonzero nuclear spin J, the quadratic Stark shift of the state \"/JFMp) is 



AW E (7, J, F, Mp , E) = -Ia scalar ( 7 , J, F)E 2 - ±a tcnsor ( 7 , J, F) 



[3M 2 - F(F + 1)] 
F(2F- 1) 



(3£ 2 2 - E 2 



(19) 



We make the approximation that hyperfine interac- differential Stark shift between the hyperfine levels of 
tion does not modify the electronic part of the atomic the ground state, which is significant for the 199 Hg + 
wavef unctions (the JJ-coupling approximation of Angel microwave frequency standard JHJ , requires going to a 
and Sandars [[19|). This approximation is adequate for higher order of perturbation theory 20 1. In the li- 
the present purpose, which is to evaluate the Stark shift coupling approximation p9| , 
of the 199 Hg + optical clock transition. Obtaining the 



a S calar(7i Ji F) = a sca i ar (7, J), 



«tc 



r ( 7 ,J,F) = (-l) 



I+J+F 



F{2F - 1){2F + 1)(2J + 3) (2 J + 1)( J + 1)1 1/2 



(2F + 3)(F+ 1)J(2J- 1) 



F J I 
J F 2 



ate 



Equations (|l7j) and (|18|) were used to evaluate the po- 
larizabilities for the Hg+ 5d 10 6s 2 Si/ 2 and 5d 9 6s 2 2 D 5 / 2 
levels. For the calculation of a sca i ar (S, 1/2), the os- 
cillator strengths for all electric-dipole transitions con- 
necting the 5d 10 6s configuration to the 5d 10 np (n = 
6,7,8) and 5d 9 6s6p configurations were included. These 
were taken from the theoretical work of Brage et al. 

The final result is a sca i ar (S, l/2)/(47re ) = 2.41 x 
10~ 24 cm 3 , which compares very well with the value of 
2.22 x 1CP 24 cm 3 obtained by Henderson et al. from 
a combination of experimental and calculated oscilla- 
tor strengths |2^j . For the calculations of a sca i ar (D, 5/2) 
and atensor(D, 5/2), the oscillator strengths for electric- 
dipole transitions to the 5d 10 np (n = 6,7,8), 5d 10 5f, 
and 5d 9 6s6p configurations were taken from Brage et 
al. pi]] . The oscillator strengths for electric-dipole 
transitions to the 5d 9 6s7p and 5d 8 6s 2 6p configurations 
were taken from the Cowan-code calculations. The re- 
sults were a sca i ar (D, 5/2)/(47re ) = 3.77 x 1CT 24 cm 3 and 
aten S or(D,5/2)/(47re ) = -0.263 x 10" 24 cm 3 . Evalu- 
ating Eq. @ for F=2 and F=3 in the 5d 9 6s 2 2 D 5/2 



(20) 

r (7,A (21) 
r (D,5/2) and 



level, we obtain a t cnsor(D, 5/2, 2) = ^a te 

atensor 

(D,5/2,3) (D,5/2). ' 

The tensor polarizability is much smaller than the 
scalar polarizabilities and in any case does not contribute 
if the external electric field is isotropic, as is the case for 
the blackbody radiation field. The net shift of the op- 
tical clock transition due to the scalar polarizabilities is 
5[a sca iar(S, 1/2) - o sca i ar (D, 5/2)].E 2 . In frequency units, 
the shift is — 1.14 X 10~ 3 E 2 Hz, where E is expressed 
in V/cm. The error in the coefficient is difficult to esti- 
mate, particularly since it is a difference of two quanti- 
ties of about the same size. However, the total shifts are 
small for typical experimental conditions. If the electric 
field is time-dependent, as for the blackbody field, the 
mean-squared value (E 2 ) is taken. At a temperature of 
300 K, the shift of the optical clock transition due to the 
blackbody electric field is —0.079 Hz. The mean-squared 
blackbody field is proportional to the fourth power of 
the temperature. For a single, laser-cooled ion in a Paul 
trap, the mean-squared trapping electric fields can be 
made small enough that the Stark shifts are not likely to 
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be observable 



5. ELECTRIC QUADRUPOLE SHIFT 

The atomic quadrupole moment is due to a depar- 
ture of the electronic charge distribution of an atom from 
spherical symmetry. Atomic quadrupole moments were 
first measured by the shift in energy levels due to an 
applied electric-field gradient in atomic-beam resonance 
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experiments 

The interaction of the atomic quadrupole moment 
with external electric-field gradients, for example those 
generated by the electrodes of an ion trap, is analo- 
gous to the interaction of a nuclear quadrupole moment 
with the electric field gradients due to the atomic elec- 
trons. Hence, we can adapt the treatment used for the 
electric-quadrupole hyperfine interaction of an atom [p6| . 
The Hamiltonian describing the interaction of external 
electric-field gradients with the atomic quadrupole mo- 
ment is 



H = V£ (2) ■ ( 



2 

E 

q=-2 



{-\yVE^Q% (22) 



where VE^ is a tensor describing the gradients of the 
external electric field at the position of the atom, and 
is the electric-quadrupole operator for the atom. 
Following Ref. J2(| , we define the components of V£' 2 ^ 

as 

1 9£ 2 

(23) 



(2) 



2 dz 



= ±f ^± = ±f d± E z 
6 az 6 



VE. 



(2) 
±2 



V6 
12 



d±E±, 



(24) 
(25) 



where E± = E x zb \E y and d± 



i).r 



(2) 

The operator components <d q arc defined in terms of 



the electronic coordinate operators as 
e 
2 



Bo P0 =-fE( 3 ^-^ 



(2) _ 
U ±l — 



e 



(2) 
±2 




( Xj ±iyj) 2 , 



(26) 



(27) 



(28) 



where the sums are taken over all the electrons. The 
quadrupole moment 0("f, J) of an atomic level \jJ) is 
defined by the diagonal matrix element in the state with 
maximum Mj: 

&( 7 ,j) = { 7 jj\e^\ 1 jj). (29) 

This is the definition used by Angel et al. [^4| . 

In order to simplify the form of V-E*- 2 "*, we make a 
principal-axis transformation as in Ref. |27[]. That is, we 



express the electric potential in the neighborhood of the 
atom as 

y', z') = A[{x' 2 + y> 2 - 2z' 2 ) + e(x' 2 - y' 2 )}. (30) 

The principal-axis (primed) frame (x',y',z') is the one 
in which <I> has the simple form of Eq. (p0|), while the 
laboratory (unprimed) frame (x,y,z) is the in which the 
magnetic field is oriented along the z axis. 

The tensor components of V£' 2 ' in the principal-axis 
frame are obtained by taking derivatives of &(x' , y' , z'): 



(2)' 



-2A 



VE (2 l' = 0, 



v4f - 



eA. 



(31) 
(32) 

(33) 



In the principal-axis frame, Hq has the simple form 



H 



Q 



-2AQ 



(2)' 



-de 



(34) 



As long as the energy shifts due to Hq are small rel- 
ative to the Zeeman shifts, which is the usual case in 
practice, Hq can be treated as a perturbation. In that 
case, it is necessary only to evaluate the matrix elements 
of Hq that are diagonal in the basis of states \-fJFMp), 
where F is the total atomic angular momentum, includ- 
ing nuclear spin J, and Mp is the eigenvalue of F z with 
respect to the laboratory (not principal-axis) frame. Let 
lj denote the set of Euler angles {a,/3, 7} that takes the 
principal-axis frame to the laboratory frame. To be ex- 
plicit, starting from the principal-axis frame, we rotate 
the coordinate system about the z axis by a, then about 
the new y axis by j3, and then about the new z axis by 7 
so that the rotated coordinate system coincides with the 
laboratory coordinate system. We can set 7 = 0, since 
the final rotation about the laboratory z axis, which is 
parallel to B, has no effect. The states \ jJFm)' defined 
in the principal-axis frame and the states defined 
in the laboratory frame are related by 



\7JFm)' = YDW( U )\-yJFri 



(35) 



( F) 

where D^ni(u;) is a rotation matrix element defined in 

E9|. The inverse relation 



the passive representation 
is 



\ 1 JF^)=J2 D i F r,rn\lJFmY. 



(36) 



In order to evaluate the diagonal matrix elements of 
Hq in the laboratory frame, it is necessary to evaluate 

(2) 1 

matrix elements of the operators <d q , defined in the 
principal-axis frame. These matrix elements are of the 
form 
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{ 1 JF^\Qf\ 1 JF^)= ]T d£,MD$>) '< 7 J^m'|e( 2 )'|7^m)', 

m' m 

= { 1 jF\\e^\\ 1 jF)Y J (-i) 



\F-r 



= (-1) 



( 7 jF||e( 2 )|| 7 jF) £ 



(-l) F -^«(7JF||G( 2 )|| 7 JF) ^ 

K m m' n n 



F 2 F 
—m q m 

F 2 F 
-m! q m 

(2K+1) 



<^)<T("), 
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,(<*>)> 
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I — m 




m / 






n' J 


\ m' 


—m 


n J 



(37) 
(38) 
(39) 



i?S*( W ),(40) 



= (-l)^-«( 7 JF||e( 2 )|| 7 JF)^^ I JjDfl/H, 



(41) 

), and @ follow from Eqs. (4.2.7U4.3.2), 

1 (with 



where Eq. (|38|) follows from the Wigner-Eckart theorem, and Eqs. (|39|), 

and (3.7.8) of Ref. [p8| , respectively. The required rotation matrix elements are, from Eq. (4.1.25) of Ref 
correction of a typographical error), 

£>$*(w) = i(3cos 2 /3-l), 



D, 



'00 

(2) 
±2 



The 3-j symbol in Eq. rt4J) is 



F 2 F 
—fj, /i 



= (-1) 



F-u 



sin 2 /3(cos 2a =F i sin 2a) . 



2[3// 2 -F(F + l)] 



[(2F + 3)(2F + 2)(2F + 1)2F(2F - 1)]V2 ' 
The diagonal matrix elements of i?Q in the laboratory frame are 
( 7 JFMp | i?Q 1 7 JFMp ) = 



(42) 
(43) 

(44) 



-2[3M 2 - F(F + l)]A( 7 JF||e< 2 )|| 7 JF) 
[(2F + 3)(2F + 2)(2F + 1)2F(2F - I)] 1 / 2 

It is simple to show, by directly integrating the angu- 
lar factor in square brackets in Eq. (|45|) , that the average 
value of the diagonal matrix elements of Hq, taken over 
all possible orientations of the laboratory frame with re- 
spect to the principal-axis frame, is zero. This also fol- 
lows directly from the fact that the quantity in square 
brackets is a linear combination of spherical harmonics. 
It is less obvious that the average, taken over any three 
mutually perpendicular orientations of the laboratory z 
quantization axis, is also zero. This result is proven in the 
Appendix. This provides a method for eliminating the 



[3 cos 2 /3 - 1) - e sin 2 /3(cos 2 a - sin 2 a) 



(45) 



quadrupole shift from the observed transition frequency. 
The magnetic field must be oriented in three mutually 
perpendicular directions with respect to the trap elec- 
trodes, which are the source of the external quadrupole 
field, but with the same magnitude of the magnetic field. 
The average of the transition frequencies taken under 
these three conditions does not contain the quadrupole 
shift. 

The reduced matrix element in Eq. ( ^5| ) is, in the IJ- 
coupling approximation, 



( 7 (F/)F||e( 2 )|| 7 (F/)F) = (-1)'+^(2F+1) 





2 




( J 2 




{i 






^ — J 


i) 



e( 7 ,j), 



(46) 



where / is included in the state notation in order to spec- 
ify the order of coupling of / and J. For the particular 
case of the 199 Hg + 5d 9 6s 2 2 D 5 / 2 level, the reduced matrix 
elements are 



(D5/2 2||e (2) ||D5/2 2) = 2 J — 6(D,5/2), (47) 

V 5 



(D5/2 3||e (2) ||D5/2 3) = 2 J— 9(D,5/2). (48) 

V 5 

Since the Cowan-code calculation shows that there is 
very little configuration mixing in the 199 Hg + 5<i 9 6s 2 
2 D 5 / 2 level, 0(D,5/2) can be reduced to a matrix ele- 
ment involving only the 5d orbital: 
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0(0,5/2) = - (5d 2 d 5/2 , mj = 5/2|3z 2 -r 2 |5d 2 d 5/2 , = 5/2), (49) 

= | (5d, mi = 2|3z 2 -r 2 |5d, m ; = 2), (50) 

= ey — (5d, m, = 2|r 2O (0,0)|5d, m ; = 2), (51) 

= e yy(5d|r 2 |5d) y y F 2 y^)y 2 o(M)r 2 2(M)sinM^, (52) 



= 5e(5d|H|5d)(^_ 2 o sJ^q o J, (53) 

= -y(5d|r 2 |5d). (54) 

The apparent sign reversal in Eq. ( fl9] ) relative to Eqs. ( |26| ) and ( p9| ) is due to the fact that the quadrupole moment 
is due to a single hole in the otherwise filled 5c? shell rather than to a single electron. According to the Cowan-code 
calculation, 

(5d\r 2 \5d) = 2.324 a 2 Q = 6.509 x 10~ 17 cm 2 , (55) 

where ao is the Bohr radius. 

Since the quadrupole shifts are zero in the 5d 10 6s 2 S]y 2 level, the quadrupole shift of the 199 Hg + optical clock 
transition is due entirely to the shift of the |D 5/2 2 0) state, and is given by 

(D 5/2 2 0|if Q |D 5/2 2 0) = f A6(D, 5/2)[(3 cos 2 j3 - 1) - esin 2 /3(cos 2 a - sin 2 a)}, (56) 

= -^Ae(5d|r 2 |5d)[(3cos 2 /3- 1) - e sin 2 /3(cos 2 a - sin 2 a)], (57) 
w -3.6 x 10- 3 hA[(3 cos 2 /3 - 1) - e sin 2 /3(cos 2 a - sin 2 a)] Hz, (58) 

where A is expressed in units of V/cm 2 . Thus, for typical values A « 10 3 V/cm 2 and |e| < 1, the quadrupole shift is 
on the order of 1 Hz. 
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6. APPENDIX. ANGULAR AVERAGING OF THE QUADRUPOLE SHIFT 

For the purpose of describing the quadrupole shift, the orientation of the laboratory (quantization) axis with respect 
to the principal-axis frame is defined by the angles (3 and a. In the principal-axis coordinate system, a unit vector 
along the laboratory z axis is defined in terms of (3 and a by 

z = (sin [3 cos a, sin (3 sin a, cos (3) . (59) 

We wish to show that the angular dependence of the quadrupole shift is such that the diagonal matrix elements given 
by Eq. ( |45"| ) average to zero, for z along any three mutually perpendicular directions. 

An arbitrary set of three mutually perpendicular unit vectors e\, e 2 , and e% can be parameterized by the set of 
angles 9, <fi, and ip in the following way: 

e\ = (sin(9cos0,sin#sin</>, cos 6), (60) 
e 2 = (cos (f>cos 9 cos ip — sin 4> sin ip, sin <j> cos 9 cos ip + cos cf> sin ip, — sin 9 cos ip), (61) 
e3 = (— cos cos 9 sin-0 — sin cos ip, — sin 4> cos 9 sin-0 + cos (/>cos?/>, sin 9 sin -0). (62) 

It can be verified by direct computation that ej • e,j = Sij . 
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The quadrupole shift can be evaluated for each of these three unit vectors substituted for z [Eq. (g^)] and the 
average taken. First consider the average of the quantity (3 cos 2 ,9—1) that appears in Eq. ( filf ) : We use the fact that 
cos (3 is the third component of z, so the average is: 



(3 cos 2 (3-1) = cos 2 9 + sin 2 9 cos 2 t/> + sin 2 9 sin 2 — 
= cos 2 9 + sin 2 6-1, 
= 0, 



(63) 
(64) 
(65) 



for arbitrary 9, cf>, and ip. Similarly, the average of the other angle-dependent term in Eq. (f45|), sin 2 /3(cos 2 a — sin 2 a), 
is calculated by making use of the fact that sin /3 cos a is the first component of z, and sin f3 sin a is the second: 



1 [sin 2 9 cos 2 ( 



sin 2 9 sin 2 i 



(sin /3(cos a — sin a)) — ., 

+(cos 4> cos 9 cos — sin sin i/>) 2 — (sin </> cos cos ip + cos </> sin ip) 2 
+ (cos</> cos 9 sin V> + sin</>cos'0) 2 — (sin cos 6 1 sin ip — cos (p cos ip) 2 } 
= 0, 



(66) 
(67) 



for arbitrary 9, <f>, and ip. Hence, the matrix elements of 
Hq given by Eq. ( fi5| ) average to zero for any three mu- 
tually perpendicular orientations of the laboratory quan- 
tization axis. 
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